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PROCEDE DE RECONSTRUCTION D 1 UNE IMAGE RADIOGRAPH I QUE 
PAR COMBINAISON DE VIGNETTES SE RECOUVRANT 

DESCRIPTION 

5 

L* invention concerne un procede de 
reconstruction d'une image radiographique par 
combinaison d'une collection de vignettes se 
recouvrant . 

10 Les images de grande taille en radiographie 

sont freguemment obtenues par morceaux, au moyen de 
vignettes de projection qui sont prises successivement 
en deplagant le rayonnement traversant l'objet a 
differentes positions, ainsi que le reseau de 

15 detecteurs prenant les mesures.. L' assemblage des 
vignettes donne ensuite 1 ' image recherchee . 

Le procede se complique dans le cas usuel 
d ! un rayonnement divergeant d'un foyer vers le reseau 
de detecteurs, soit en cone, soit en un ensemble 

20 d'eventails plans et paralleles. La figure 1 montre la 
configuration normale des mesures : le rayonnement 
comprend une source 1 (ponctuelle ou lineaire) qu 1 on 
deplace a chaque mesure le long de l'objet 2 ainsi que 
le reseau de detecteurs 3. Les positions prises sont 

25 notees la, lb, lc, Id, et 3a, 3b, 3c et 3d. Pour que 
1 ' attenuation du rayonnement soit mesuree a tout 
endroit de l'objet 2, le faisceau 4 du rayonnement doit 
comprendre des portions de recouvrement dans les 
positions 4a, 4b, 4c et 4d qu 1 on lui fait prendre, et 

3 0 les vignettes de projection du rayonnement, dont les 
positions coincident avec celles 3a, 3b, 3c et 3d que 
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le reseau de detecteurs 3 prend successivement , 
presentent de meme des portions de recouvrement . II est 
done impossible de juxtaposer simplement les vignettes 
pour obtenir 1 ' image globale de l'objet, mais on doit 
au contraire determiner les positions des portions de 
recouvrement sur les vignettes et faire une synthase du 
contenu de ces portions de recouvrement pour 
reconstruire 1 1 image . 

Un autre probleme qui apparait est celui du 
grossissement des details selon leur distance a la 
source 1. La largeur de projection de details 5 de 
l'objet 2 sur le reseau de detecteurs 3 sera 
proportionnellement plus large si les details 5 sont 
plus proches de la source 1. Un rayonnement divergent 
ne permet ainsi ni de juxtaposer facilement des 
vignettes, ni de respecter l'echelle des details a 
l'interieur de chaque vignette. 

La figure 2 permet d'expliciter ces 
problemes. Deux details 5a et 5b espaces verticalement 
se trouvent dans la portion de recouvrement des 
vignettes prises par le reseau de detecteur 3 aux 
positions 3a et 3b. Les rayons passant par le detail 5a 
sont distants de 1 ' ecart 6 sur le plan du reseau de 
detecteurs 3, et ceux qui traversent le detail 5b sont 
distants de 1 1 ecart 7 sur le meme plan ; les ecarts 
sont differents entre eux, et differents du deplacement 
qu'il a fallu imposer au reseau de detecteurs 3 entre 
les positions 3a et 3b ou les vues ont ete prises. Une 
bonne reconstruction de 1 ' image aux portions de 
recouvrement impose de combiner entre elles les mesures 
associees a chacun des details pour des vignettes 
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dif f erentes , ce qui est impossible a faire directement 
puisque leurs hauteurs sont inconnues en general . Si on 
choisit par exemple d'associer les rayons distants de 
1'ecart 6 pour reconstruire les portions de 
5 recouvrement, les details a hauteur de 5a seront rendus 
correctement, mais les details presents a d'autres 
hauteurs ne pourront pas I'etre. La combinaison des 
vignettes produira alors du flou et un grossissement 
inexact de ces autres details . 

10 Un procede de ce genre a cependant deja ete 

propose dans 1 ' art . L' image est reconstruite en 
choisissant une reconstruction exacte des portions de 
recouvrement a des hauteurs determinees, la ou les 
details important s, et notamment les os pour une 

15 radiographic du corps, sont susceptibles de se trouver. 
Pour obtenir des bons resultats, il faut done connaitre 
d'abord la hauteur de ces details. Des coefficients de 
ponderation peuvent privilegier les resultats de Pune 
ou 1' autre des vignettes selon la position consideree 

20 sur la portion de recouvrement. La restitution des 
autres details de 1 ' image est sacrifiee. 

Un autre procede connu consiste a calculer 
des correlations entre les portions de recouvrement des 
dif f erentes vignettes pour evaluer 1'ecart (6, 7 ou 

25 autre) des rayons a associer pour faire la synthese des 
portions de recouvrement. Les correlations dependent de 
details preponderant s presents sur les deux portions de 
recouvrement et provenant d'un meme endroit de l'objet 
2. La reconstruction de 1 ' image est accomplie a hauteur 

3 0 de ces details preponderants et elle est bonne, si du 
moins ces details existent ; mais comme dans le procede 
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precedent, les details situes aux autres hauteurs 
seront mal rendus . 

II faut a j outer que des conflits de hauteur 
peuvent apparaitre si les portions de recouvrement sont 
norabreuses, et en particulier avec un rayonnement 
conique ou les portions de recouvrement concernent tout 
le perimetre des vignettes. Deux portions de 
recouvrement sur deux cotes d'une vignette pourront 
etre reconstruites independamment a des hauteurs 
differentes, tout en ayant une intersection pour 
laquelle on sera embarrasse a choisir une hauteur de 
reconstruction. 

Un procede de reconstruction plus correcte 
d'une image radiographique est propose avec 
1' invention. II est fonde sur une discretisation 
generale de l'objet en volumes (voxels) definissant des 
hauteurs de reconstruction, et des combinaisons des 
valeurs d' attenuation estimees sur chacun des volumes 
aux differentes hauteurs de reconstruction pour 
ameliorer 1' image globale, sans privilegier 

necessairement une hauteur de reconstruction. 

De facon plus detaillee, 1' invention 
concerne generalement un procede de reconstruction 
d'une image radiographique d'un objet traverse par un 
rayonnement divergent subissant une attenuation, le 
rayonnement occupant des positions successives ayant 
des portions de chevauchement et 1 • attenuation etant 
mesuree par un reseau de detecteurs, sur lesquels le 
rayonnement se projette et donnant des vignettes de 
1- image respect ivement associees aux positions du 
rayonnement et comprenant aussi des portions de 
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chevauchement , le procede comprenant une combinaison de 
vignettes pour reconstruire 1 ' image , ainsi que les 
etapes suivantes : 

discretiser l'objet en volumes 
5 definissant des hauteurs de reconstruction, 

associer les volumes a au moins un 
detecteur respectif du reseau sur lequel le rayonnement 
se projette apres avoir traverse ledit volume, 

attribuer une valeur d ■ attenuation a 

10 chaque volume d 1 apres les valeurs mesurees par ledit 
detecteur associe , 

- et combiner les valeurs d 1 attenuation des 
volumes aux differentes hauteurs de reconstruction. 

Dans une des formes de 1' invention, la 

15 valeur d' attenuation attribuee a chaque volume est 
egale a la somme des valeurs mesurees par - ledit 
detecteur associe, divisee par le nombre des vignettes 
que contribuent a donner ledit detecteur associe, et 
les valeurs d r attenuation des volumes sont combinees 

2 0 par une combinaison numerique sur des groupes des 
volumes superposes aux differentes hauteurs de 
reconstruction. Et dans une autre de ses formes, la 
valeur d' attenuation attribuee a chaque volume est 
obtenue par retroproj ect ion iterative des valeurs 

25 d 1 attenuation mesurees par les detecteurs, des valeurs 
provisoires etant attributes aux volumes et corrigees 
apres avoir ete projetees sur les detecteurs, en 
calculant des differences entre des sommes et les 
valeurs provisoires sur des lignes de projection aux 

30 valeurs mesurees par les detecteurs sur lesdites lignes 
de projection, et en repartissant les differences sur 
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lesdites lignes de projection pour corriger les valeurs 
provisoires . 

L' invention sera raaintenant decrite 
completement en liaison aux figures, dont la figure 1 
schematise le precede de realisation des vignettes, la 
figure 2 illustre le probleme de reconstruction a une 
hauteur arbitraire, la figure 3 illustre les elements 
participant a 1 ' exploitation de 1' invention, et les 
figures 4 et 5 sent des organigrammes de deux modes du 
procede . 

On passe au commentaire de la figure 3 . 
L'objet 2 est discretise en volumes ou voxels 8 qui 
definissent des hauteurs de reconstruction 11. L e 
rayonnement passe a travers les volumes 8 par des 
rayons 9i et 9 j , qui sont plusieurs pour les volumes 8 
appartenant aux portions de recouvrement , et qui sont 
originaires de positions respectives li et lj de la 
source 1 et se projettent sur des detecteurs respect if s 
10i et lOj qui leur sont associes pour les positions 
correspondantes 3i et 3j du reseau de detecteurs 3. Les 
detecteurs 10 mesurent des attenuations des rayons 9 a 
travers tout l'objet 2, et done a travers tous les 
volumes tels que 8 qu'ils traversent . Dans la pratique, 
les volumes 8 se projettent sur une surface qui peut 
englober plusieurs detecteurs 10 completement, et 
d'autres partiellement . Le systeme est calibre pour 
associer a chaque volume 8 sa surface de projection et 
calculer 1 ' attenuation qu'il fait subir au rayonnement 
en proposition des detecteurs concernes . Nous 
n'evoquerons pas ici ces techniques de calibration, qui 
sont tout a fait usuelles dans l'art, et considererons 
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des volumes 8 se projetant completement sur un 
detecteur 10 unique le long d'un rayon de projection 
unique par souci de simplicity des explications. 

D'apres la figure 4, le procede commence 
5 done par une etape A de discretisation generale de 
1'objet 2 en volumes 8 dont les couches definissent les 
hauteurs de reconstruction 11 de 1 ' image . En pratique, 
les hauteurs de reconstruction 11 seront assez peu 
nombreuses et les volumes 8 plutot parallelepipediques , 

10 allonges en hauteur, que cubiques. Les etapes suivantes 
B et C consistent a se placer a une hauteur de 
reconstruction 11 et un volume 8 determines. On 
recherche alors les rayons tels que 9i et 9j passant 
par le volume 8 considere, et les detecteurs lOi et lOj 

15 de projection desdits rayons sur le reseau de 
detecteurs 3, a 1' etape D. L' etape suivante E consiste 
a lire la mesure d' attenuation des rayons 9i et 9j sur 
les detecteurs lOi et 10 j . A 1' etape suivante F, une 
moyenne de ces attenuations est faite, du raoins pour 

2 0 les volumes 8 appartenant aux portions de recouvrement 

et qui sont done traverses par au moins deux rayons 9 . 
Les etapes C a E ou F sont ensuite repetees pour tous 
les volumes de la couche consideree ; apres quoi, a 
1' etape G, une image de 1'objet 2 est reconstrui te . 
25 Cette image est une image de 1' ensemble de 

1'objet 2, et pas seulement une image de coupe a la 
hauteur consideree, puisque les attenuations mesurees 
par les detecteurs 10 le long des rayons 9 ont ete 
supposees concentrees aux volumes 8 de la couche a 

3 0 cette hauteur. 
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Ensuite, on revient a l'etape B pour 

reconstruire l'objet 2 a une autre hauteur, et le cycle 

des etapes C a G recommence avec les volumes 8 de la 
couche associee. 

Quand les images de l'objet 2 ont ete 
reconstruites a toutes les hauteurs, elles sont 
combinees a l'etape H avec l'espoir d'obtenir une image 
plus exacte. Plusieurs precedes peuvent etre envisages. 
Le plus simple consiste peut-etre a faire des moyennes 
des images sur des colonnes 12 (a la figure 3) de 
volumes 8 empiles appartenant a des couches 
differentes, avec eventuellement une ponderation pour 
favoriser les couches les plus representatives. 
Eventuellement, on peut choisir une seule des images 
qu'on juge meilleure que les autres, ou un assemblage 
de plusieurs des images aux endroits qu' elles 
representee le mieux. Toutes ces methodes devraient 
donner de meilleurs resultats que celles de 1 ' art 
anterieur qu'on a decrites auparavant . 

On ne fera que mentionner certains procedes 
de correction qui sont usuels dans 1 ' art et qui ne sont 
pas affectes par 1' invention. 

Le rayonnement diffuse peut tout d'abord 
etre soustrait des mesures avant d' exploiter celles-ci. 
Plusieurs procedes existent pour faire cette 
soustraction, dont le plus simple est peut-etre de 
proceder a une mesure supplemental^ ou un ecran est 
intercale entre l'objet 2 et le reseau de detecteurs 3 
en masquant certains des detecteurs 10. Les detecteurs 
10 masques ne sont pas touches par le rayonnement 
direct des rayons 9, mais seulement par le rayonnement 
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diffuse, qui est alors mesure par ces detecteurs et qui 
peut etre deduit par des interpolations pour les autres 
detecteurs. 

Les attenuations d'un rayonnement peuvent 
5 en general etre exprimees par un coefficient 
multiplicateur du rayonnement initial Io inferieur a 
1' unite et egal a e -ix±, ou 1 est la longueur 
d' attenuation et [l le coefficient d' attenuation 
caracteristique du materiau, et qui est en general la 

10 valeur qu'on recherche pour reconstruire 1 ' image „ Les 
detecteurs 10 mesurent directement le rayonnement I qui 
n'a pas ete absorbe par 1'objet 2 et qui est egal a 
I 0 e^ 1 ; on peut en deduire le produit Jill , puis la 
valeur de |a si on divise les valeurs de |ll par les 

15 longueurs de traversee de 1'objet 2 par les rayons 9, 
apres les avoir estimees par une autre mesure ou les 
avoir evaluees. 

Un autre de mode de realisation de 
1' invention sera maintenant decrit au moyen de la 

20 figure 5. Apres une etape de discretisation J semblable 
a celle A de la realisation precedente, on effectue au 
mieux une division en blocs a 1 ' etape K. En effet, la 
resolution qui va etre entreprise peut devenir malaisee 
si le systeme considere est trop volumineux. En 

2 5 pratique chaque bloc pourra comprendre les volumes 8 
associes a une vignette. Qu'une division en bloc soit 
faite ou non, le probleme a resoudre peut s'exprimer 
par p = Mx ou x designe les inconnues, c'est-a-dire 
les attenuations aux volumes 8, p designe les 

30 projections de ces valeurs, c'est-a-dire les mesures 
par les detecteurs 10, et enfin M designe la matrice de 
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projection. Les coefficients de la matrice M 

represented la contribution d'un volume 8 d' indice j a 
la projection suivant le rayon 9 d' indice i, et peuvent 
en general etre approches par la longueur traversee par 
ce rayon dans ce volume. 

L'etape suivante L est une evaluation de 
1' attenuation aux volumes 8 du bloc considere. La 
premiere evaluation pent etre arbitraire, par exemple a 
des valeurs nulles. Pour chacun des volumes 8 on 
recherche, d'apres l'etape M, le detecteur 10 qui lui 
est associe par le rayon 9 qui l e traverse, comme a 
l'etape D de la realisation precedente. L'etape 
suivante N est une lecture des mesures des detecteurs 
10 semblable a l'etape E. La determination des rayons 9 
de projection permet de proceder a une evaluation des 
valeurs projetees d' attenuation a l'etape O, c'est-a- 
dire qu'on precede au calcul Mx pour evaluer p. En 
soustrayant ces valeurs evaluees des projections aux 
valeurs reelles, mesurees, des memes projections, on 
determine 1'erreur commise dans 1' evaluation des 
valeurs projetees a l'etape P. 

L'etape suivante Q est une retroprojection 
de cette erreur dans les volumes 8 de 1'objet 2 afin de 
corriger les valeurs evaluees de 1 ' attenuation 
Concretement, on precede en executant la formule 

INw^df bloc " M ""* q ) 00 * <q+ " ^ S « sont 
des evaluations successives de 1 ' attenuation aux 
volumes 8 du bloc ; X<« est un coefficient de 
relaxation permettant de ne pas se diriger trop vite 
vers une solution qui ne correspond qu'aux premiers 
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blocs et qui est compris entre 0 et 2 ; ce coefficient 
n'est d'ailleurs pas uniforme dans les blocs mais peut 
avantageusement etre plus eleve pour les rayons 
sensiblement verticaux, ou perpendiculaires aux 
5 detecteurs 10, afin de leur donner une importance plus 
grande ; t Mbioc est la transposee de la matrice M pour le 
bloc considere ; le terme au denominateur est un terrae 
de normalisation ; enfin, les termes entre parentheses 
representent 1 ' erreur calculee a l'etape P . 

10 On procede de meme pour le bloc suivant, en 

recommengant le cycle de l'etape K a l'etape Q, puis on 
revient au premier bloc pour une nouvelle iteration, 
jusqu'a ce que les attenuations evaluees aient converge 
vers une solution, ce qu'on exprime par l'etape R . Les 

15 volumes 8 inclus dans les portions de recouvrement du 
volume 2 ont ete traites de la meme fagon que les 
autres, en subissant simplement des iterations plus 
nombreuses qu'eux s'ils sont associes a plusieurs blocs 
de projections. 

2 0 On dispose alors d'une image 

tridimensionnelle de l'objet 2 ; une image 
bidimensionnelle de bonne qualite peut etre obtenue par 
une combinaison des valeurs obtenues, qui consiste a 
additionner les valeurs d' attenuation sur les colonnes 

25 12 de volumes 8 empiles. 

Le procede de 1' invention permet de 
concilier une bonne qualite de restitution des details 
importants de l'objet etudie avec une bonne qualite 
d' ensemble de 1 ' image . II est possible d'obtenir des 

30 images dont la resolution est analogue au pas des 
detecteurs 10. 
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On s'est place dans la situation usuelle ou 
le reseau 3 de detecteurs accompagne le mouvement du 
rayonnement 4, mais le precede pourrait etre applique 
sans changement avec un reseau de detecteurs immobile 
sous 1'objet 2 et dont la surface s'etendrait a toutes 
les vignettes de projection. 



3 14215o3/JCI 



1er depot 

13 

REVEND I CAT IONS 

1) Procede de reconstruction d'une image 
radiographique d'un objet traverse par un rayonnement 
5 divergent subissant une attenuation, le rayonnement 
occupant des positions successives (4) ayant des 
portions de chevauchement et 1 1 attenuation etant 
mesuree par un reseau (3) de detecteurs (10), sur 
lesquels le rayonnement se projette et donnant des 

10 vignettes de 1 1 image respect ivement associees aux 
positions du rayonnement et comprenant aussi des 
portions de chevauchement, le procede comprenant une 
combinaison de vignettes pour reconstruire l f image, 
ainsi que les etapes suivantes : 

15 - discretiser 1 'objet en volumes (8) 

def inissant des hauteurs de reconstruction (11), 

associer les volumes a au moins un 
detecteur respectif du reseau sur lequel le rayonnement 
se projette apres avoir traverse ledit volume, 

2 0 - attribuer une valeur d 1 attenuation a 

chaque volume d ! apres les valeurs mesurees par ledit 
detecteur associe , 

- et combiner les valeurs d ' attenuation des 
volumes aux dif ferentes hauteurs de reconstruction. 
25 2) Procede de reconstruction d 1 une image 

radiographique selon la revendicat ion 1, caracterise en 
ce que la valeur d T attenuation attribuee a chaque 
volume est egale a la somme des valeurs mesurees par 
ledit detecteur associe, divisee par le norabre des 

3 0 vignettes que contribuent a donner ledit detecteur 

associe, et les valeurs d 1 attenuation des volumes sont 
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combinees par une combinaison numerique sur des groupes 
(12) des volumes superposes aux differentes hauteurs de 
reconstruction . 

3) Procede de reconstruction d'une image 
radiographique selon la revendication 1, caracterise en 
ce que la valeur d 1 attenuation attribute a chaque 
volume est obtenue par retropro j ection iterative des 
valeurs d 1 attenuation mesurees par les detecteurs (10) , 
des valeurs provisoires etant attributes aux volumes et 
corrigees apres avoir ete projetees sur les detecteurs, 
en calculant des differences entre des sommes des 
valeurs provisoires sur des lignes de projection et les 
valeurs mesurees par les detecteurs sur lesdites lignes 
de projection, et en repartissant les differences sur 
lesdites lignes de projection pour corriger les valeurs 
provisoires . 

4) Procede de reconstruction d'une image 
radiographique selon la revendication 3, caracterise en 
ce que les valeurs d 1 attenuation des volumes sont 
combinees en une somme sur des groupes (12) des volumes 
superposes aux differentes hauteurs de reconstruction. 



25 
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